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ABSTRACT 

We investigate the properties of multiple-image gravitational lens systems formed by 
dark matter halos. Both the Navarro-Frenk- White (NFW) and Moore et al. mass profiles 
are considered, and the effects of quadrupole perturbations to the lensing potential are 
taken into account. The systems produced by dark halos exhibit two generic properties 
that serve as powerful benchmarks for scrutinizing individual lens candidates: small 
flux ratios between the two brightest images in triples, and detectable odd images. In 
contrast, most currently ambiguous quasar pairs consist of two components with flux 
ratios > 3:1. Such systems are statistical outliers in the expected distribution of dark 
lenses, and are therefore likely to be binary quasars. 

Subject headings: gravitational lensing - dark matter - quasars 

1. Introduction 

Observational evidence indicates that virtually all galaxies are embedded in dark matter halos, 
but do all halos host galaxies? Differences in the slope of the mass function for galaxies (e.g., Blanton 
et al. 2001) and dark halos (Press & Schechter 1974) suggest that a substantial fraction of halos 
may be empty, particularly at lower masses (e.g., Nagamine et al. 2001). Because cold dark matter 
(CDM) simulations predict cuspy halos (e.g., Navarro, Frenk & White 1997; Moore et al. 1999b; 
Jing & Suto 2000) that are sufficiently concentrated to multiply image background sources (Li &, 
Ostriker 2002; Wyithe, Turner &: Spergel 2001), strong gravitational lensing may be a powerful 
tool for detecting empty dark matter halos. 

The identification of gravitational lens candidates is based on the presence of two or more 
closely-separated quasars at the same redshift. In cases where a galaxy, group or cluster is respon- 
sible for the multiple imaging, the lens hypothesis can be unambiguously confirmed by directly 
detecting the lensing mass through optical or infrared observations. This is clearly not possible 
for dark lenses. A thorough investigation of the quasar components is therefore necessary to de- 
termine whether a system is a genuine dark lens, or a pair of physically distinct quasars. The 
lensing hypothesis can be robustly rejected by finding some spectral or photometric property that 
the quasars do not share in common. For example, a pair with only one radio-bright component 
(an 2 R system, see Kochanek, Falco & Munoz 1999) is inconsistent with gravitational lensing, 
and must therefore be a binary quasar. The inability to uncover differences in the quasars is not 
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compelling evidence for lensing, however, as similarities can be at least partially attributed to the 
common formation history and cosmic environment of binaries. Unless the presence of three or 
more quasar components makes the lensing identification obvious, a measurement of the time delay 
or a detection of correlated milliarcsecond-scale radio substructure would likely be necessary to 
prove that a system is a dark lens. 

A number of wide-separation (AO > 3") quasar pairs with no detectable lensing mass have 
been reported (e.g., Hawkins 1997; Mortlock, Webster & Francis 1999), and their true nature has 
been the subject of continuing debate. Several of the quasar pairs are 2 R systems (e.g., MGC 
2214+3550; Munoz et al. 1998), and therefore binaries. For others the lensing hypothesis has been 
rejected based on incompatible spectral line widths (e.g., CTQ 839; Morgan et al. 2000) or X-ray 
properties (e.g., Q2345+007; Green et al. 2002). Yet, there exist ~ 10 systems for which the 
identification remains ambiguous. Statistical arguments suggest that most of these pairs are likely 
to be binary quasars (Kochanek et al. 1999). The scarcity of radio-bright /radio-bright dark lens 
candidates is particularly difficult to reconcile with a sample dominated by lenses, considering how 
many radio lenses have been discovered. Still, the existence of several dark lenses is allowed by 
statistical constraints, and such arguments cannot determine the nature of any one source. It is 
therefore desirable to find additional criteria by which to scrutinize individual dark lens candidates. 

Most studies of strong gravitational lensing by dark matter halos have focussed on their contri- 
bution to the optical depth (Keeton & Madau 2001; Li & Ostriker 2002; Sarbu, Rusin & Ma 2001), 
rather than the properties of the lens systems produced. This paper bridges the gap, and will serve 
as a 'field guide' for dark lenses. In §2 we investigate lenses formed by plausible models of dark ha- 
los, quantifying the distribution of magnification ratios between the two brightest images in triples, 
as well as the relative flux of the third image. In §3 current dark lens candidates are compared to 
our model predictions. Section 4 summarizes our findings and discusses their implications. 

2. Predictions for Dark Lens Models 

Numerical simulations predict CDM halos with cuspy but shallow inner mass profiles (p oc 
r~ a at small radii, where 1.0 < a < 1.5; Navarro et al. 1997; Moore et al. 1999b). This result 
applies to halos in their "primordial" state. Cooling baryons will modify the dark matter profile 
through adiabatic compression (e.g., Blumenthal et al. 1986; Keeton 2001; Kochanek & White 
2001), producing the more centrally-concentrated mass distributions that are characteristic of lens 
galaxies (e.g., Cohn et al. 2001; Kochanek 1995a; Rusin & Ma 2001). Inversely, if the baryons 
have not cooled, the dark matter halo is expected to retain its primordial profile (e.g., Navarro, 
Frenk & White 1995; Eke, Navarro & Frenk 1998). Star formation should be greatly inhibited in 
such halos, so the mass distributions will be truly "dark" (except perhaps at X-ray wavelengths). 
The existence of a characteristic mass scale separating cooled halos with nearly isothermal profiles 
(efficient lenses) and non-cooled halos with primordial dark matter profiles (inefficient lenses; e.g., 
Wyithe et al. 2001) is vital for converting the halo mass function into the observed distribution of 
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lensed image separations (Porciani & Madau 2000; Kochanek & White 2001). Specifically, halos 
above this scale must have a much lower cross section per unit mass to account for the sharp decline 
in observed lenses with A6 > 3". The necessary mass scale (M c ~ 1O 13 M ) agrees well with the 
predictions of semi-analytic cooling models (e.g., Cole et al. 2000). To produce arcsecond-scale 
(1" < A9 < 10") lens systems, halos that retain their shallow primordial profiles require virial 
masses M v ; r > 1O 135 M (e.g., Li & Ostriker 2002) - larger than the cooling scale. Consequently, 
the formation of observable multiple-image lenses by uncooled and largely unmodified CDM halos 
is self-consistent. Note, however, that the lack of star formation may not imply an uncompressed 
halo. In particular, baryons in halos with M v - ir < 1O 12 M can cool into a disk that is stable to 
star formation if the spin parameter is high (Jimenez et al. 1997). Such halos would be both dark 
and concentrated. But even if the resulting mass profile were as steep as isothermal, the limit of 
1O 12 M corresponds to a lensed image separation of less than an arcsecond (e.g., Li & Ostriker 
2002). The above arguments therefore suggest that arcsecond-scale dark lens systems are likely to 
be produced by shallow CDM mass profiles. 

This paper investigates two popular models for cuspy CDM halos: the Navarro-Frenk- White 
(NFW; 1997) mass distribution (a = 1), and the steeper (a = 1.5) profile of Moore et al. (1999b). 
While the inner profile slope may be mass dependent (Jing & Suto 2000), the above models span the 
typical range of values found. Spherical mass models are considered for computational simplicity, 
but we place them in constant shear fields (magnitude 7) to mimic the effects of deflector ellipticity 
and environmental perturbations. We focus on two/three- image configurations, which we term 
"triples," because these are the systems most likely to be confused with binary quasars. The 
deflector models are analyzed by numerically solving for the image positions and magnifications 
(Hi) of point sources placed on a uniform grid behind the lens. The distributions of two quantities 
are then extracted: the magnification ratio between the two brightest images (r±2 = \ni/n2\), and 
that between the brightest and third brightest images (ri3 = \p,i/fis\). 

The properties of dark lenses differ significantly from those of the more concentrated, nearly 
isothermal (p oc r~ 2 ) galaxy lenses. Specifically, the magnification ratio distributions are very 
sensitive to the size of the radial critical curve, which increases as the slope of the inner mass 
profile is decreased from isothermal (Blandford & Kochanek 1987; Wallington & Narayan 1993; 
see Kormann, Schneider & Bartelmann 1994 for plots). First, source positions enclosed by only 
the radial caustic produce triples in which the second (negative-parity) image resides outside of the 
radial critical curve; the central (positive-parity) image resides within it. Large radial critical curves 
imply that the outer images tend to form at similar distances from the lens center, resulting in 
similar magnifications and small primary flux ratios r±2 . The third image can also form farther from 
the lens center, away from the high central convergence that would otherwise strongly demagnify 
it. Consequently, this image carries more flux, leading to smaller values of 7-13. Second, the critical 
structures of shallower mass profiles are more susceptible to shear (e.g., Wallington & Narayan 
1993), which allows for the production of naked cusp configurations. Here, source positions enclosed 
by only the tangential caustic are lensed into three images on the same side of the deflector, each 
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with similar magnification (e.g., Kormann et al. 1994). This class of triples is characterized by very 
small magnification ratios and bright third images, thereby reinforcing the above trends. 

Flux ratio distributions are affected by magnification bias, which favors the detection of more 
magnified lens systems (e.g., Turner 1980). The bias factor depends on the lens model, the number- 
flux relation of the sources, and the method of finding lenses. A full accounting of magnification bias 
is beyond the scope of this paper. However, we illustrate the bias effect by simply weighting each lens 
system generated on the source plane grid by its total magnification (/j, = J2i | Mi I ) ■ For example, 
the biased fraction of lenses with r 12 < r is [/ y 2 )R(yi, y 2 , r 12 )dy 1 dy 2 ]/[J n(y\,y 2 )dy\dy 2 \, 
where R = 1 if r\ 2 < vq for a source at (yi,y 2 ), and otherwise. The integral is evaluated over the 
three-image caustic regions. This scheme mimics the bias produced by the differential number-flux 
relation N(S) oc S~ 2 , which closely approximates the population of compact radio sources (Rusin 
& Tegmark 2001). Faint optical quasars have a much shallower luminosity function, and hence 
they are nearly unbiased (Kochanek 1996). Bright optical quasars are more strongly biased, as the 
luminosity function is very steep in that regime. 



2.1. NFW Profile 

The NFW model is described by the mass density 

, x_ Pcr^I (0 -.s 

P(r) ~ {r/R s )[l + {r/R s )Y ^ 

where R s = R 2 oo/c(z), R 2 oo is the radius within which the average density is 200 times the critical 
density p cr at redshift z, and c(z) is the concentration parameter. For a flat 17a = 0.7 cosmology, 
co = c(0) ~ 7 (e.g., Bartelmann et al. 1998) and evolves as c(z) = cq/(1 + z) (e.g., Bullock et al. 
2001). The density amplitude is <5j = (200/3)c 3 /[m(l+c)-c/(l+c)]. The spherical NFW profile has 
an analytic projected surface density and deflection angle (Bartelmann 1996). Its lensing properties 
depend on the parameter kq = pcrSiRgfEcr, where S cr is the critical surface mass density. The 
angular image splitting scale (A6>) is set by the size of the tangential critical curve. In Fig. 1, A9 
is plotted as a function of kq for some typical redshifts and concentrations. 

The fractions of NFW triples with a magnification ratio between the two brightest images of 
ri2 < 1.5, 2.0 and 5.0 are plotted in Fig. 2. For clarity, we plot the cumulative distributions of r\ 2 
for kq = 0.10, 0.15, and 0.20 in Fig. 3. Note that magnification bias only slightly alters the values, 
and does not affect the general trends. For a spherically-symmetric potential (7 = 0), almost all 
triples have r\ 2 ~ 1. Shear perturbations create a tangential caustic on the source plane, affecting 
the distribution of magnification ratios. If the cusp of the tangential caustic resides within the radial 
caustic, a single highly magnified image results for sources near the cusp. Such systems can have 
large r\ 2 , and therefore the fraction of lenses with r\ 2 less than some nominal value is decreased. 
If the cusp is not enclosed by the radial caustic, however, the deflector can produce naked cusp 
configurations. Such systems have r\ 2 < 2, and therefore the fraction of lenses with small r\ 2 is 
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increased. Deflectors with smaller kq are more susceptible to shear perturbations, as the ratio of 
the radial to tangential critical curve increases with decreasing kq (Bartelmann 1996). Hence, for 
a fixed 7, the fraction of systems with small 772 rises sharply below that hq at which naked cusp 
configurations begin to dominate. This effect is clearly seen in Figs. 2 and 3. The conclusion is 
that NFW dark lenses with 1" < A6 < 10" (using the conversion in Fig. 1) are expected to have 
flux ratios very close to unity if either the potentials are nearly spherical, or have moderate to 
large shear perturbations (7 > 0.10). But even in the intermediate regime, systems with high flux 
ratios are still rare. For 7 = 0.05, nearly half of arcsecond-scale NFW lenses will have 772 < 2, and 
almost all will have 772 < 5 (Fig. 2). In comparison, of the lenses produced by a singular isothermal 
sphere, 11% have 772 < 2 and 44% have 772 < 5. 

The relative magnification of the third image (773 = I/J1//X3Q is examined in Figs. 4 and 5. 
Plotted in Fig. 4 are the fractions of NFW triples with 773 < 20, 50 and 100. Plotted in Fig. 5 
are the cumulative distributions of 773 for kq = 0.10, 0.15, and 0.20. The third image tends to be 
very bright and easily detectable, independent of kq or 7. In almost all cases, at least half of the 
triples will have 773 < 20. The third image is even more detectable in regions of parameter space 
dominated by naked cusps, as these systems have three images of very similar magnification. 



2.2. Moore Profile 

The model of Moore et al. (1999b) is described by the mass density 

/ \ = Per S II . . 

P{r) (r/R 9 )W[l + (r/R a )W] { ~> 

where Su = 100c 3 / ln(l + c 3 / 2 ), and c is a factor of a few smaller than the NFW value. There is 
no analytic form for the projected surface density or deflection angle. Typical input parameters, 
however, produce lensed images at radii that are much smaller than the arcminute-scale angular 
break radius derived from R s . Because lenses will be dominated by the inner cusp, the Moore profile 
can be well represented by the power-law mass density p(r) oc r _L5 . The critical structures of power- 
law models scale in a self-similar manner with the profile normalization, and hence the lensing 
properties are normalization- independent. Trial analyses of the real Moore profile demonstrate 
that it is well approximated by the power law and exhibits the expected normalization behavior. 
We therefore present results for the singular power-law model, which were calculated using the 
deflection angles and magnification matrices derived by Barkana (1998) and implemented in the 
FASTELL software package. 

Plotted in Fig. 6 are the cumulative distributions of 772 and 773 for the Moore profile. The 
unbiased fractions of Moore triples with 772 < 1.5, 2.0 and 5.0 are, respectively, 0.99, 1.00 and 1.00 
for 7 = 0; 0.45, 0.71 and 0.97 for 7 = 0.05; 0.44, 0.69 and 0.98 for 7 = 0.10; and 0.59, 0.77 and 0.95 
for 7 = 0.15. The unbiased fractions of Moore triples with 773 < 20, 50 and 100 are, respectively, 
0.32, 0.49 and 0.61 for 7 = 0; 0.33, 0.51 and 0.64 for 7 = 0.05; 0.44, 0.69 and 0.86 for 7 = 0.10, and 
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0.74, 0.99 and 1.00 for 7 = 0.15. The largest difference between the Moore and NFW predictions 
for arcsecond-scale triples is due to the naked cusp configurations. The tangential caustic is less 
prominent for a Moore model than it is for an NFW model in the same shear field. Consequently, 
while naked cusps comprise a significant fraction of NFW triples for small shears (7 < 0.05), they 
do not contribute at all to the Moore statistics until the shear grows large (7 > 0.10). But despite 
this effect, dark halos described by the steeper mass profile of Moore et al. still produce lensed 
triples with small magnification ratios and relatively bright third images. 

3. Comparison to Current Candidates 

A number of dark lens candidates have been reported in the literature, and their true nature 
is still a matter of vigorous debate (e.g., Hawkins 1997; Kochanek et al. 1999; Mortlock et al. 1999; 
Peng et al. 1999). Table 1 lists nine ambiguous quasar pairs. We have excluded those systems for 
which the lensing hypothesis has been rejected based on gross spectral or photometric differences 
between the components. Technically Q2345+007 (Weedman et al. 1982) is no longer ambiguous, 
as recent X-ray observations have shown it to be a binary (Green et al. 2002), but we include it 
for historical reasons. Below we determine whether the component fluxes of these quasar pairs are 
compatible with the dark lens hypothesis. 

First, dark lens candidates typically consist of two, not three, detectable images. Bounds on 
additional images in these systems are listed in Table 1. The absence of third images down to these 
levels (ri3 5 ii m it > 50) would be very unlikely in a candidate sample dominated by lenses, particularly 
if dark halos were best approximated by an NFW model (Figs. 4 and 5). Second, the observed flux 
ratios in many of the ambiguous pairs are significantly larger than unity. Three of the nine pairs 
have ri2,obs > 10, which Figs. 2, 3 and 6 show to be virtually impossible for genuine dark lenses. 
However, even the quasar pairs with moderate flux ratios (3 < f72 i0 bs < 5) would be uncommon 
in the dark lens hypothesis. Note that the paucity of small flux ratio pairs cannot be due to the 
selection function, as their identification is always favored (Kochanek 1995b). 

Quantitative results are presented in Table 2. For each quasar pair and model, we calculate the 
unbiased fraction of triples that have a primary magnification ratio 772 smaller than the observed 
value, or a third image brighter than the listed bound. Fractions with the third image constraint 
excluded are also shown. Each of the systems is an outlier in the predicted dark lens distributions, 
particularly those for the NFW profile. Five of the nine pairs have p(r\2 < r72, bs V?73 < r^rm^) > 
0.95 in all eight models studied. These pairs are strongly incompatible with dark lensing. Three 
other pairs have 0.90 < p(r\2 < ?72,obs V 773 < 773^^) < 0.95 in each of the models. QJ0240-343, 
with an observed flux ratio of 772 = 2.1, is moderately compatible with one of the four Moore models, 
but no NFW models. Magnification bias has little effect on the above results. In conclusion, while 
no quasar pair can be absolutely rejected as a dark lens based on probability arguments alone, the 
above analysis demonstrates that most current candidates fail to exhibit the properties expected 
for dark lenses. Therefore, these pairs are likely to be binary quasars, as suggested by Kochanek 
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et al. (1999). 

The only system that meets the qualitative dark lens criteria is APM 08279+5255 (Ibata et al. 
1999), which consists of three quasar components with a maximum separation of 0?.'38. The outer 
two components carry nearly equal flux (ri2,obs — 1-3), and the central component is very bright 
( r i3,obs — 5.7). The striking similarity of the quasar colors and the presence of a third compact 
component means that APM 08279+5255 is almost certainly a lens. No lensing galaxy has yet 
been detected, but this is not a surprise considering the brightness of the quasars (m# ~ 13 — 14). 
Radio observations have recently detected line emission offset from the position of the three lensed 
components (Lewis et al. 2002), which suggests imaging by a naked cusp. Because third images are 
rare for the steep mass profiles of galaxy lenses (Wallington & Narayan 1993; Rusin & Ma 2001), 
the unique morphology of APM 08279+5255 opens the possibility of strong lensing by a shallow 
dark matter halo. Assuming that the system is a standard triple (source inside the radial caustic), 
Muhoz, Kochanek & Keeton (2001) have demonstrated that a cuspy profile even shallower than 
NFW is necessary to reproduce the image data. Our Monte Carlo tests find that the flux ratios 
of APM 08279+5255 are moderately consistent with an NFW model normalized to a present-day 
concentration of cq ~ 7— 10, but only if the quadrupole is small (7 < 0.05). Naked cusps dominate 
for larger shears, and such systems tend to have values of r±2 and ri3 much smaller than the 
observed ratios. In contrast, while Moore triples commonly form with 1.3 < r\i < 2.0 (Fig. 6), 
the third image tends to be significantly fainter. Hence, while APM 08279+5255 exhibits the 
qualitative properties expected for a dark lens, none of the standard models provide a particularly 
good quantitative match to the system. 



4. Discussion 

This paper has investigated the properties of gravitational lens systems produced by dark 
matter halos. While galaxy lenses form with a wide range of magnification ratios and typically 
have undetectable central images, dark lenses are expected to have very small magnification ratios 
and prominent third images. These characteristics result from the shallow inner mass profiles of 
dark halos, as derived by N-body simulations. Relative component fluxes in candidate dark lens 
systems are powerful benchmarks for evaluating their viability. The majority of current candidates 
have flux ratios that differ significantly from unity, and do not feature any third image. These 
sources therefore fail to exhibit the expected properties of genuine dark lenses. While the above 
arguments cannot absolutely reject any source as a dark lens, they offer further evidence that binary 
quasars dominate the sample. The only system that meets the qualitative dark lens criteria is APM 
08279+5255 (Ibata et al. 1999). 

There are several possible criticisms of the results presented in this paper. First, the dark 
matter profiles were assumed to be smooth, despite the fact that simulations predict copious sub- 
structure in CDM halos (e.g., Moore et al. 1999a). Substructure can alter the magnifications of 
individual lensed images (Mao & Schneider 1998), an effect that may have already been detected 
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(Dalal Sz Kochanek 2002). However, Monte Carlo tests indicate that substructure tends to decrease 
magnification ratios in doubles (Metcalf & Madau 2001), at least for isothermal profiles. If this 
result can be extrapolated to the shallow halo profiles, then high flux ratio dark lenses may be 
even rarer than we predict. Second, we assumed that dark lenses will have concentration parame- 
ters similar to those of the general halo population. The simulations of Bullock et al. (2001), for 
example, predict a range of halo concentrations, and lensing will tend to select more concentrated 
mass distributions because of their larger cross sections. Will this make the gravitational potentials 
of dark lenses significantly steeper? Consider NFW halos with concentration parameters a factor 
of three larger than their typical values. There are not many halos found in this range (Bullock 
et al. 2001). Lenses with 1" < A9 < 10" then correspond to normalizations of 0.25 < kq < 0.55 
(0.3 < Zd < 0.7). Because Figs. 2 and 4 show that such lenses still have very small primary flux ra- 
tios and bright third images, "concentration bias" should not alter our general conclusions. Finally, 
one might suggest that the measured flux ratios may not be good estimators of the true magnifi- 
cation ratios if these pairs were real lenses, because of either differential extinction, microlensing, 
or a conspiracy between source variability and time delays. This is unlikely to be a major concern, 
however. Similarity between the quasar colors (e.g., Tyson et al. 1986; Hewett et al. 1998), often 
cited in support of lensing, argues against the extinction hypothesis. Furthermore, the relatively 
small levels of variation observed in the flux ratios - in some cases over a baseline of more than a 
decade - suggest that time delays and microlensing are not significantly clouding the measurements. 

While we have focussed on triples, five-image lens systems are expected to play an important 
role in dark lensing. Small quadrupole perturbations will result in a large tangential caustic for 
shallow profiles, thereby producing many quads. Plotted in Fig. 7 are the unbiased fractions of 
five-image lenses formed by the NFW model. More than half of NFW lenses may be quads, even for 
very small shears. These fractions will be further increased by magnification bias. Using the Moore 
profile, the unbiased quad fractions are 0.06, 0.25 and 0.47 for 7 = 0.05, 0.10 and 0.15, respectively. 
The absence of quads makes the current dark lens candidate sample even more peculiar. 

Dark lenses would be powerful probes of the mass distributions in primordial dark matter halos 
- assuming they can be found. Unfortunately, the shallow profiles mean that CDM halos are very 
inefficient lenses (e.g., Wyithe et al. 2001). Cases of strong lensing are thus expected to be rare, 
even if empty halos were abundant. Consequently, the best hope for uncovering dark lenses are 
large systematic quasar surveys. None were found in the Cosmic Lens All-Sky Survey (e.g., Myers 
et al. 1999; 2002), as lens galaxies have been identified in virtually all of the systems (Jackson et al. 
1998; Browne et al. 2002), and radio spectral data strongly suggest that the quasar pair B0827+525 
(Koopmans et al. 2000) is a physical binary. The Sloan Digital Sky Survey, however, offers new 
hope in the search. The expected properties of dark lenses outlined in this paper should be a useful 
tool for critiquing any candidates identified therein. 

The author is grateful to Chris Kochanek, Chung-Pei Ma, Nick Sarbu and Stuart Wyithe 
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Table 1. Historically Ambiguous Quasar Pairs 



Source 


-^comp 


Type 


Zs 


Ad 


^12,obs 


n3,limit 


Reference 


MG0023+171 


2 




0.95 


A'!8 


~ 10 a 


~ 20 


Hewitt et al. 1987 


QJ0240-343 


2 


O 2 


1.41 


&!l 


2.1 


~ 100 6 


Tinney 1995 


B0827+525 


2 


2 R 2 


2.06 


2?8 


2.8 C 


~ 100 


Koopmans et al. 2000 


Q1120+0195 


2 


O 2 


1.47 


6'.'5 


> 60 


> 100 


Meylan & Djorgovski 1989 


LBQS1429-008 


2 


O 2 


2.08 


5'/l 


17.4 


> 100 


Hewett et al. 1989 


Q1634+267 


2 


O 2 


1.96 


3"8 


3.3 d 


~ 50 


Djorgovski & Spinrad 1984 


Q2138-431 


2 


O 2 


1.64 


4'/5 


3.0 


~ 70 


Hawkins et al. 1997 


LBQS2153-2056 


2 


O 2 


1.85 


T!8 


14.5 


> 100 


Hewett et al. 1998 


Q2345+007 


2 


O 2 


2.15 


7'!3 


4.0 


> 100 e 


Weedman et al. 1982 



Note. — Listed for each pair are the number of quasar components (N comp ), radio-brightness 
classification, source redshift (z s ), angular separation (A9), observed flux ratio (r 12j0 bs)) bound on 
additional components (ri3 limit ) 5 an d reference. Many of the parameters are taken from the tables 
of Peng et al. (1999) or Mortlock et al. (1999). Radio flux ratios are less susceptible to microlensing 
and extinction, and are used where available. Values of ri3 } n m i t are from the listed references, except 
where noted. a 5 GHz measurement (optical ratio is ~ 3); b Cf A- Arizona Space Telescope Lens 
Survey; c 8.4 GHz measurement; d updated value from Peng et al. (1999); e Tyson et al. (1986). 
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Table 2. Quasar Pairs vs. the Dark Lens Distribution 









NFW profile 










Moore profile 






Source 


7 = 


0.05 


7 = 


0.10 


7 = 


0.15 


7 = 


0.05 


7 = 


0.10 


7 = 


0.15 


MG0023+171 


0.99 


(0.99) 


1.00 


(0.99) 


1.00 


(1.00) 


1.00 


(1.00) 


1.00 


(1.00) 


0.99 


(0.99) 


QJ0240-343 


1.00 


(0.63) 


1.00 


(0.65) 


1.00 


(1.00) 


0.85 


(0.74) 


0.96 


(0.72) 


1.00 


(0.78) 


B0827+525 


1.00 


(0.76) 


1.00 


(0.92) 


1.00 


(1.00) 


0.91 


(0.85) 


0.97 


(0.86) 


1.00 


(0.86) 


Q1120+0195 


1.00 


(1.00) 


1.00 


(1.00) 


1.00 


(1.00) 


1.00 


(1.00) 


1.00 


(1.00) 


1.00 


(1.00) 


LBQS 1429-008 


1.00 


(1.00) 


1.00 


(1.00) 


1.00 


(1.00) 


1.00 


(1.00) 


1.00 


(1.00) 


1.00 


(1.00) 


Q1634+267 


1.00 


(0.83) 


1.00 


(0.92) 


1.00 


(1.00) 


0.92 


(0.90) 


0.96 


(0.93) 


0.99 


(0.90) 


Q2138-431 


1.00 


(0.80) 


1.00 


(0.87) 


1.00 


(1.00) 


0.91 


(0.87) 


0.96 


(0.90) 


0.99 


(0.88) 


LBQS2153-2056 


1.00 


(1.00) 


1.00 


(0.99) 


1.00 


(1.00) 


1.00 


(1.00) 


1.00 


(1.00) 


1.00 


(1.00) 


Q2345+007 


1.00 


(0.90) 


1.00 


(0.91) 


1.00 


(1.00) 


0.95 


(0.94) 


0.98 


(0.96) 


1.00 


(0.93) 



Note. ■ - Listed for each quasar pair and model are the predicted fractions of triples p(r±2 < 
fi2,obs V ri3 < r^iimit) and, in parentheses, p(r\2 < ?"i2,obs)- These probabilities are 1.00 (1.00) 
for all pairs in spherical (7 = 0) models. For the NFW profile, the fractions were calculated by 
deriving kq from the angular splitting scale A9, assuming a lens redshift z,i = 0.5, a source redshift 
z s = 2.0 (measured redshifts were ignored for simplicity), a present-day concentration parameter 
Co = 7, and a flat S7a = 0.7 cosmology. Varying these parameters within their reasonable ranges 
has a negligible effect on the results. All the model predictions are unbiased. Including bias has 
little qualitative impact on the derived fractions. 
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Fig. 1. — Conversion between the NFW normalization (kq) and angular splitting scale (A0) for 
fixed source redshift z s = 2, present-day concentration parameters cq = 7 (solid) and 10 (dash), 
and lens redshifts Zd = 0.3, 0.5 and 0.7 (left to right). A flat S7a = 0.7 cosmology is assumed. The 
separation A9 is approximated as twice the size of the tangential critical curve for the spherical 
case. Naked cusp configurations tend to have somewhat smaller image separations than standard 
configurations produced by the same deflector. However, the small correction does not affect the 
general trends discussed in this paper, and is therefore ignored. Note that arcsecond-scale lenses 
are dominated by 0.1 < kq < 0.2. 
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7 = 



y = 0.05 unbiased 



P(r- 
P(r- 
p( r - 



< 1.5) 

< 2.0) 

< 5.0) 



7 = 




- 7 - 0.10 


unbiased 












y = 0.05 


biased 










y = 0.15 biased 




Fig. 2. — Primary magnification ratio ryi = I/X1//X2I in NFW triples: fractions (p) of triples with 
7*12 < 1.5 (solid), 2.0 (dash) and 5.0 (dot-dash) as a function of kq for shear amplitudes 7 = 0, 0.05, 
0.10 and 0.15. Unbiased and biased values are shown. 




Fig. 3. — Distribution of magnification ratios r\i in NFW triples: fractions p{r\2 < ?"o) f° r K o = 0.10 
(solid), kq = 0.15 (dash) and kq = 0.20 (dot-dash). Unbiased and biased values are shown for shear 
amplitudes 7 = 0, 0.05, 0.10 and 0.15. 
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v d 



7 = 



p(r 13 < 20) 
p(r 13 < 50) 
p(r 13 < 100) 



7 = 



7 = 0.05 unbiased 




7 = 0.10 unbiased 



7 = 0.15 unbiased 



7 = 0.05 biased 



7 = 0.10 biased 



7 = 0.15 biased 



Fig. 4. — Same as in Fig. 2, but for the brightness of the third image in NFW triples (ri3 = |/ii/^3| 
fractions (p) of triples with r\s < 20 (solid), 50 (dash) and 100 (dot-dash). 



-7 = 


unbiased 




k = 0.10 

k = 0.15 

k = 0.20 



- y = 0.05 


unbiased 


/ " 
/ '' 




/ '' 
















Fig. 5. — Same as in Fig. 3, but for ratio r±s: fractions p(r±s < ro). 
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Fig. 6. — Distribution of magnification ratios r±2 (left) and r\z (right) in Moore triples. Plotted are 
the fractions p(ri2,i3 < tq) for 7 = (solid), 0.05 (dash), 0.10 (dot-dash) and 0.15 (dot). Unbiased 
and biased values are shown. 
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Fig. 7. — Fraction of five-image systems produced by the NFW model as a function of kq for 
7 = 0.05 (dash), 0.10 (dot-dash) and 0.15 (dot). Only unbiased values are shown. Note that the 
curves turn over for small kq. In this regime the number of quads is diluted, as some sources within 
the tangential caustic are now lensed into cusp configurations. Consequently, deflectors with kq to 
the left of the above peaks produce triples that are dominated by naked cusps. 



